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We introduce a minimal theory of glass formation based on the ideas of molecular crowding and 
resultant string-like cooperative rearrangement, and address the effects of free interfaces. In the bulk 
case, we obtain a scaling expression for the number of particles taking part in cooperative strings, and 
we recover the Adam-Gibbs description of glassy dynamics. Then, by including thermal dilatation, 
the Vogel-Fulcher-Tammann relation is derived. Moreover, the random and string-like characters 
of the cooperative rearrangement allow us to predict a temperature-dependent expression for the 
cooperative length ^ of bulk relaxation. Finally, we explore the influence of sample boundaries 
when the system size becomes comparable to The theory is in agreement with measurements 
of the glass-transition temperature of thin polymer films, and allows to quantify the temperature- 
dependent thickness hm of the interfacial mobile layer. 


According to Philip Anderson, the deepest and 
most interesting unsolved problem in solid-state 
physics is probably the glass transition. By exten¬ 
sion, this includes the highly debated confinement 
effects in glassy polymer films. The present ar¬ 
ticle introduces a minimal analytical model, that 
invokes only the ideas of molecular crowding and 
string-like cooperative rearrangement, before ad¬ 
dressing the key effects of interfaces. The valid¬ 
ity and simplicity of the approach make it ideal 
for application to various systems and geometries, 
and suggest that dynamics in glass-forming ma¬ 
terials might be understood from elementary ar¬ 
guments. 

Glassy materials are ubiquitous in nature [T], and dis¬ 
cussions about glass transition involve many areas of 
physics, from molecular and spin glasses to hard-sphere 
jamming PH?]. In spite of the intense interest in the 
dynamical slowing that accompanies glass formation, a 
single microscopic theory has yet to emerge [8HT5]. Nev¬ 
ertheless, the phenomenological approach of free vol¬ 
ume [14] and the Doolittle ansatz m have been used to 
support the Vogel-Fulcher-Tammann (VFT) relation [iBl - 
im, which describes so many of the observed behav¬ 
iors. Fundamental to glass formation are the suggestions 
that particles are increasingly crowded, and relaxation 
requires the cooperative participation of a growing num¬ 
ber of particles. The hypothesis of a cooperatively rear¬ 
ranging region, as introduced by Adam and Gibbs [H] , is 
appealing and has been observed in computational stud¬ 
ies pniiiij- 

The existence of a length scale ^ for cooperative re¬ 
arrangement |22j has led to tremendous interest in con¬ 
fined glass formers, as initiated by [23 • Perhaps, the 
most active example of attempts to probe ^ is the study 
of glassy polymer films P^26j . where fascinating obser¬ 


vations have been made. For the most studied case of 
polystyrene, reductions in the measured glass-transition 
temperature have been almost uniformly reported as the 
film thickness is reduced, both experimentally mi and 
numerically [28]. It has been further suggested that this 
apparent anomaly is linked to the observed existence of a 
more mobile interfacial layer P51l32j . As a consequence, 
there have been many theoretical attempts to understand 
the thin-film glass transition, with varying degrees of 
complexity and success [33H38] . 

In this article, we present a simple analytical model 
for relaxation in glass-forming materials. First, from a 
microscopic molecular picture, the nature of the cooper¬ 
ative mechanism is explicitly defined and characterized as 
a function of density, and the Adam-Gibbs phenomenol¬ 
ogy is recovered. Then, by including thermal expansivity, 
we derive the VFT relation for the temperature depen¬ 
dence of the relaxation time in bulk materials. Finally, 
in order to address the effects of interfaces, the theory is 
applied to the case of thin films. 

Beyond any formulation, there are two main ingredi¬ 
ents that a microscopic cooperative theory must contain: 
i) “more cooperative is easier”, and ii) “more cooper¬ 
ative is rarer”. The first one means that, in order to 
redistribute a given amount of volume in a crowded en¬ 
vironment, a cooperative rearrangement is energetically 
more favourable than a solitary one - because the for¬ 
mer is the sum of N small displacements which is easier 
to satisfy than a single large displacement. This effect 
tends to maximize N. The second ingredient relies on 
the fact that glass-forming materials are made of inde¬ 
pendent particles that move incoherently due to thermal 
fluctuations. Therefore, it is relatively rare to have mo¬ 
tions that are coherent in time and space and form col¬ 
lective objects: the larger N, the rarer the event. This 
effect tends to minimize N. Taken together, these two 
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FIG. 1: Two relaxation modes in glass-forming materials, (a) 
At low volume fraction ij), or high temperature T, a test par¬ 
ticle (green) can escape its neighbouring cage by single parti¬ 
cle motion, (b) At higher volume fraction, or lower tempera¬ 
ture, the previous mechanism is inhibited, but relaxation of the 
test particle can still occur through a random string-like co¬ 
operative process: neighbours (blue) that temporarily get into 
close proximity with each other provide additional space for 
the relaxation. At the end of the string, an incoherent parti¬ 
cle (grey) terminates the process. 


ingredients suggest a most probable value of N. Chang¬ 
ing the temperature may change the crowding constraints 
and/or the coherence penalty, which result in the temper¬ 
ature dependence of this most probable value and thus 
the glassy behavior of interest. From numerical simula¬ 
tions [201 HU HD comes another important feature: it has 
been observed that the cooperative regions often have a 
fractal dimension close to 1. This was also reported in 
experimental studies of repulsive colloids [30]. Therefore, 
we add a third ingredient to the list above: iii) “cooper¬ 
ative rearrangement is string-like”. In the following, our 
goal is to build the simplest mean-field toy-model that 
contains i), ii), and iii). 

As shown in Fig. we consider an assembly of parti¬ 
cles of effective radius r, average intermolecular distance 
A, and volume fraction cj) ex (r/A)^. A given particle is 
surrounded by a cage of nearest neighbours. Relaxation 
requires passage between two adjacent neighbours, the 
“gate”, with average size L ^ X — 2r. We define A = Ay, 
with Ay ^ 2r, as the point of kinetic arrest with volume 
fraction (j)y. Note that, because we are only considering a 
single type of motion, and because the pressure is finite, 
the value of i/y will be less than the jamming volume 
fraction c/j m- As density increases, prior to complete 
kinetic arrest, there exists an onset value (/>c < </>v of 
the volume fraction, and associated Ac > Ay of the in¬ 
terparticle distance, at which cooperative rearrangement 
becomes the only possible relaxation mechanism. Note 
that for a cooperative motion to actually be a cooperative 
rearrangement, at least two particles have to exchange 
their positions. Note also that the value of Ac can be 


very close to Ay in actual physical systems, for several 
reasons such as deformability [41] and anisotropy [42| of 
the molecules. 

As a reference, we consider the liquid-like case of a par¬ 
ticle escaping from its cage by solitary motion (Fig. 

On average, such motion is allowed if A > Ac, by def¬ 
inition of Ac. When A > Ac, the probability density 
of relaxation per unit time and per unit volume is thus 
Pc ~ 1/(tc A^), where the constant Tc is a typical liquid¬ 
like relaxation time at the cooperative onset. At the cur¬ 
rent level of minimal description, the Boltzmann factor 
associated with the sharp repulsive intermolecular poten¬ 
tial of the gate has been replaced by an implicit Heaviside 
function on A — Ac. 

When A < Ac, solitary escape cannot occur since the 
gate is too small by an average length 5 ~ Ac — A. Re¬ 
laxation is only possible through a cooperative process 
involving N — 1 neighbours of the test particle, that get 
into close proximity with each other. Thereby, the miss¬ 
ing space 6 can be locally and temporarily made avail¬ 
able, thus allowing for a rearrangement, as observed in 
bidisperse hard disks [13]. Inspired by computational 
studies [iniHiii30] and experiments m, we consider co¬ 
operative regions in the form of string-like random chains 
(Fig#). Since the gate between particles has an average 
length L, the additional length created by the coopera¬ 
tive move is A ~ {N — 1) L. The escape of the test 
particle through a collective motion is thus possible if 
A > 6. The threshold corresponds to N = N*, where: 


iV*(</) 


Ac — Ay 
A — Ay 



( 1 ) 


This central quantity, that we will refer to as the bulk co- 
operativity, is plotted in Fig. Note that Eq. Q could 
have been assumed, independently of any details on the 
string-like microscopic picture, since it is the simplest 
expression having the expected properties: it equals 1 at 
the cooperative onset (pc', and it diverges at the kinetic 
arrest point (p\. The so-called bulk glass-transition point 
lies somewhere in between (pc and t/y, and depends 
on the time scale of the observations. Naturally, such 
a collective scenario requires all the individual motions 
making up the cooperative string to be in phase. This 
coherence is unlikely to be spontaneously satisfied with 
random isotropic molecular fluctuations, resulting in a 
penalty. The latter can be translated into the probabil¬ 
ity factor ~ (1 _ for the independent motions of 

the N — 1 consecutive particles of the cooperative string 
to occur coherently with the motion of the first test par¬ 
ticle. Here, e is the elementary coherence probability to 
be determined; and the termination factor I — e expresses 
the incoherence of the N -\- 1*^ particle (Fig. [^). Finally, 
the two features above can be combined to express the 
probability density of a cooperative relaxation process 
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FIG. 2: Bulk cooperativity as a function of normalized inverse 
volume fraction (Eg. We fixed (fv/f’c to 

realistic values for polystyrene, using a = = 5.5 x 

10-4 ^-1 j,hnik ^ 371 p. ^ 403 and 

Tv ^ 322 K (Fig. gp. 

involving N particles: 

P7v(<^) ~ ^ (1 - e) 0 (^ - N*) ^ (2) 

where, as in the previous liquid-like case, the Boltzmann 
factor associated with the sharp repulsive intermolecular 
potential has been replaced by an Heaviside function 0. 
Note that Pn{4>) is maximal for N = N*, at fixed (j). 

Summing the over all N > N*, one obtains the 
total probability density of relaxation: 

P(<^) ^ Pc (3) 

The relaxation is entirely determined by the coopera¬ 
tivity N*, and is exponentially decaying with increasing 
N*, as reported in jamming studies of air-driven gran¬ 
ular beads |44j. Introducing the ergodic correspondence 
between the bulk relaxation time r and Eq. ([^ , through 
rPA^ ~ I, and defining the molecular time scale Tq = ere, 
one obtains: 



This equation expresses the fact that the cooperative re¬ 
laxation is a combination of N* independent motions, 
that are all similar to a solitary escape at the coopera¬ 
tive onset. Note that in a thermal description, one would 
estimate the onset relaxation time through the Arrhenius 
law: Tc ~ Toe^B^c ^ where A/i is the gate energy barrier 
and k^Tc is the thermal energy at the cooperative onset. 
We thus obtain the Adam-Gibbs phenomenology Bi- 
Having described the effect of crowding on the string¬ 
like cooperativity, we now study the glass transition of 


bulk systems. In particular, we characterize the re¬ 
laxation time as a function of temperature T, by cou¬ 
pling the previous density-based picture to the ther¬ 
mal expansion coefficient a = — ^ ^ of the equilibrium 
melt state. We assume that the material dilatation is 
small in the considered range, that goes from coopera¬ 
tive onset (0c, Tc) to kinetic arrest (0v, TV)- This is 
valid for several glass formers, such as polystyrene for 
which: a = 5.5 X [H], Ty = 327 K gS], and 

Tc = 463 K [11117], so that 0(r) ~ 0v [1 -h a(Ty - T)]. 
Combining the latter with Eqs. and Q, one directly 
derives the VET relation [T^I^ . or equivalently the 
Williams-Landel-Ferry (WLF) relation [48], respectively: 


t{T) To exp 


/ A A 


\T-Ty) 

~ Tc exp 


A (Tc - T) 


(T-Ty)(Tc-Ty) 

(5) 

where Ty is identified as the Vogel temperature [TB] , and 
where A = (Ty — Tc) ln(e) is a reference temperature 
(see jH] for a thermodynamic derivation of the VET re¬ 
lation based on chain-like excitations). Therefore, the 
elementary coherence probability e = tq/tc = to/t(Tc) is 
the normalized relaxation rate at the cooperative onset. 
The bulk glass-transition temperature Tgis defined as 
the one at which r reaches a given large experimental 
time scale = r(Tg“^^). 

Since our model leads to the VET and WLF time- 
temperature superpositions, it captures well the so-called 
fragile-glass phenomenology m, and links thermal ex¬ 
pansion and fragility as observed in metallic glasses m- 
Strong-glass phenomenology can be recovered as well 
when Ty ^ Tg'^*’^. This is reminiscent of soft colloidal 
glasses [U], for which the kinetic arrest point is shifted 
to higher volume fractions due to particle deformabil- 
ity. Note that additional molecular processes, beyond 
the scope of this work, may lead to bulk relaxation that 
does not diverge at Ty |52l|53]. Besides, as the system 
approaches TV, the fractal dimension of the cooperative 
regions may change gT]. Finally, our cooperative-string 
model is a free-volume approach and, as such, is subject 
to the same criticism as all free-volume models [51ES]. 

The bulk relaxation process presented above consists 
of random cooperative strings and is entirely determined 
by N*. Therefore, in the vicinity of the kinetic arrest 
point, the length scale ^ of the cooperative regions reads 
^ ~ XN*", where: ix = 1/2 for a simple random walk; 
v « 0.588 for a self-avoiding walk; and ;/ = 1/3 for a 3D 
compact shape. Consistent with our minimal string-like 
model, we consider the simple random walk. Invoking 
Eq. 0 and thermal expansivity, one obtains: 


(iT) ~ Av . (6) 

At the cooperative onset, this asymptotic expression of 
the cooperative length corresponds to the effective par¬ 
ticle size Ay ~ 2r. At the Vogel temperature, ^ diverges 
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FIG. 3: Two string-like cooperative paths in a thin film of 
thickness h. For an inert supporting substrate (thick line), 
relaxation of a test particle (green) at a distance z from a 
free interface (dashed line) can occur through either a bulk 
cooperative string (blue) of size ( (Eq. ^), or a truncated 
string (red) touching the interface. 


with exponent —1/2. Note that this power-law is sim¬ 
ilar to the one measured in 2D vibrated granular me¬ 
dia around jamming |56j . The exponent also compares 
well to the values of —0.59 from kinetically constrained 
lattice-gas simulations m, and —0.62 from numerical 
studies of Lennard-Jones liquids [58] . Unfortunately, the 
direct measurement of ^(T) below is challenging for 
3D thermal glasses. A natural alternative way is thus to 
reduce the sample size towards ^ HSIIMIIIS]- 

We now turn to the case of a thin film of thickness h, 
supported on an inert substrate. Our model predicts that 
the key effect of the free interface is to favour a higher 
surface mobility, as observed in experiments HZiiiMsg. 
Indeed, the cooperative strings can now be truncated by 
the cage-free boundary (Fig. [^, leading to a lower sur¬ 
face cooperativity. Therefore, we introduce the average 
local cooperativity N*{z,T), at a distance z from the 
free interface. The natural length scale of the problem 
is the bulk cooperative length When z S> the free 
interface is typically not reached with less than N* co¬ 
operative particles, and the relaxation is bulk-like with 
N* ^ N*. As z —>■ 0, N* vanishes due to the absence 
of caging at the boundary. In between, N* takes all the 
intermediate values. Thus, N* varies along z over the 
bulk cooperative length and is expected to have the 
following asymptotic self-similar form in the vicinity of 
Tv: 

A/(z,r)~w, (7) 

where / is a continuous and monotonic function satisfy¬ 
ing /(O) - 0, and f{u » I) - I. 

Since our description of the cooperative process in¬ 
volves random strings of particles, Eq. Q can be sup¬ 
ported at large N* by the following argument based on 
Brownian motion. For a given cooperative string, we de¬ 
fine no < cx) as the number of particles at which the string 
hits the free interface for the first time. If uq > N*, the 
behavior is bulk-like; if no < N*, the string is truncated 


by the interface. Therefore, the quantity of interest is 
the density probability g of first passage at the inter¬ 
face with no cooperative particles. Defining the first- 
passage “time” to = no/N*, and starting at “distance” 
Z = zj^ from the interface, one gets the ID expression: 
(/(to, ^) = (27r)“^/^tQ"^^^Zexp [—Z^/(2to)] [55]. The lo¬ 
cal cooperativity being the minimum between no and 
N*, the average local cooperativity is defined as N* = 


iV* (min(l,to))t^ = N* 1 -/^ dto (1 - to) ^(to, Z) 
where we explicitly see the interfacial lowering of the 
bulk cooperativity. Calculating the integral, we recover 
Eq. Q with / (uV^) = erf(M) -I- 2uexp (—n^) /y^ — 
2u^ erfc(M). Note that the exact functional form chosen 
for / is not crucial when comparing to the experimental 
data below, as other sufficiently sharp functions provide 
similar results. 

Following the derivation of Eq. Q and assuming that 
one can replace N* by the local average cooperativity 
N* of Eq. Q, one obtains the local relaxation time: 


rs{z,T) ^ 

To [tqJ 


( 8 ) 


We thus see that / acts as a local exponent, ranging 
from 0 to 1, on the normalized bulk relaxation time. This 
formula generalizes the Adam-Gibbs phenomenology m 
by accounting for the effect of a free interface. 

Finally, we compare our theory to dilatometric mea¬ 
surements of reduced glass-transition temperatures in 
thin polystyrene films supported on silicon substrates. 
In the experiments, the thickness-dependent glass- 
transition temperature Tg{h) is defined as the location 
of a kink in the dilatation plot obtained by ellipsome- 
try |24] . This change of expansivity occurs when the 
system is half-glassy and half-liquid. Given that / is 
monotonic, this translates to the apparent transition oc¬ 
curring when, at the middle z = /i/2 of the film, the 
local relaxation time Tg equals the bulk relaxation time 
at the bulk glass transition [SD]. Invoking Eqs. (|^ 

and Tg{h) thus satisfies: 

n(r,) p) 

where f~^ denotes the inverse of the bijective function 
/. The solution of this equation is plotted in Fig. [^ 
and compared to measurements on polystyrene [611 E2] ■ 
The literature data encompasses a wide variety of tech¬ 
niques and protocols [SI], and for purposes of fitting we 
consider the restricted data of [62|, where the annealing 
conditions and atmosphere have been carefully controlled 
and documented. We see from Fig. [^ that our model pro¬ 
vides excellent agreement with the experiments. The two 
adjustable parameters are: the critical interparticle dis¬ 
tance Ay ~ 3.7 nm which is reasonably found to be com¬ 
parable to the persistent length of polystyrene [531 Hll ; 
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FIG. 4: Comparison between dilatometric experimental data 
(filled symbols) for the reduced glass-transition temperature 
Tg{h) of thin polystyrene films supported on silicon sub¬ 
strates \6^ . and the theory (line) given by Eg. Other 

literature data JSlf is shown (open symbols) for complete¬ 
ness. The fixed parameters are the bulk glass-transition tem¬ 
perature = 371 K 164'j , and the onset temperature 

Tc = 463 K fgl The two adjustable parameters are the 
critical interparticle distance Ay « 3.7 nm, and the Vogel 
temperature Ty « 322 K. 


and the Vogel temperature Ty « 322 K, which is close to 
the reported value of 327 K [35]. Note that we have only 
considered polystyrene, as that material is extremely well 
studied. However, by varying the single parameter Ay in 
Eqs. Q and (§, one would obtain larger, smaller, or even 
immeasurable reductions in 7^, as could be observed in 
other materials. 

As one notices in Fig. [^ and since f does not vary 
much around according to Eq. M, the crossover 

thickness at which the measured glass-transition temper¬ 
ature first shows deviations from the bulk value is a few 
^(j*uik)_ Interestingly, this statement is equivalent to 
a purely finite-size criterion |B5] . even though it is ob¬ 
tained from an explicit truncation of cooperative strings 
at the free interface. We may thus understand why there 
has been so much debate in the past between the purely 
finite-size effect and the mobile-layer hypothesis [55] . 

Since the comparison between theory and experiments 
(Fig.i provides an estimate of the critical interparticle 
distance Ay, we can now compute the cooperative length 
^(T) for polystyrene above and below In Fig. ^ ^ 

increases with reducing temperature and diverges at Ty, 
as more and more molecules are required to move coop¬ 
eratively for relaxation to occur. As discussed below, our 
model also predicts the existence of a liquid-like mobile 
layer at the free interface of a thin glassy film, such as 
that observed in [571 [75H57] . This introduces another rel¬ 
evant length scale to the problem: the mobile-layer thick¬ 


ness h^{T). An important matter of debate has been to 
determine the relation between ^ and hm, if any |67| . In 
other words, is a thin-film experiment able to probe the 
bulk cooperative length scale or does it introduce an¬ 
other independent length scale through the mobile-layer 
thickness h^7 An advantage of the present microscopic 
picture is to provide a tentative answer to this question. 
At a given temperature T < there exists a posi¬ 

tion z = hni where the local relaxation time Tg equals the 
bulk relaxation time at the bulk glass transition 
Invoking Eqs. ([^ and ([^, this implies: 

= . ( 10 ) 


For z > hin, Tg is larger than and the system is 

glassy; for z < h^, Tg is lower than Tg“'^ and the system 
is liquid-like. Therefore, hm is identified as the liquid-like 
mobile-layer thickness. As seen in Eq. (101 and Fig. 
hyn{T) is different but related to ^{T). Both lengths orig¬ 
inate from the same cooperative-string model, but ^ is a 
bulk quantity whereas reflects the truncation of the 
strings at the free interface. As a result, increases 
with increasing temperature and diverges at the bulk 
glass-transition temperature, when the entire material is 
in the liquid state, as observed in recent experiments m- 
To conclude, we have developed a cooperative-string 
model that connects in a predictive manner essential in¬ 
gredients of the glass transition, in bulk systems and near 
interfaces. The theory is based only on the idea of co- 
operativity required by increasing molecular crowding, 
and introduces a string-like cooperative mechanism that 
is motivated by recent studies. An outcome of our ide¬ 
alized microscopic description is to recover the Adam- 



FIG. 5: Bulk cooperative length (dotted, Eg. ^) and surface 
mobile-layer thickness hm (plain, Eq. 10)) of polystyrene, as 
a function of temperature. We used Tg = 371 K I6'4f , Tc = 
463 K I SSI \4^! the values Xv ~ 3.7 nm and Tv ~ 322 K 
obtained from the fit in Fig. 
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Gibbs picture, as well as the Vogel-Fulcher-Tammann re¬ 
lation, without the need for the Doolittle ansatz to link 
free volume and relaxation. In particular, we derive ex¬ 
plicit scaling expressions for the cooperativity and associ¬ 
ated relaxation probability. Furthermore, the simplicity 
of the model enables application to reported anomalies 
in the glass transition of thin polymer films. Specifically, 
the free interface truncates the cooperative strings and 
thus enhances the mobility in its vicinity. Agreement 
between the present theory and reported dilatometric 
measurements of the glass-transition temperature of sup¬ 
ported polystyrene films is excellent. The two adjustable 
parameters are: the critical interparticle distance at ki¬ 
netic arrest, which is found to be similar to the persis¬ 
tence length of polystyrene; and the Vogel temperature 
which is found to be close to literature values. Finally, 
the model provides a way to distinguish between purely 
finite-size and surface effects, and to clarify the existing 
link between cooperative length and mobile-layer thick¬ 
ness. Importantly, the success of the theory applied to 
the thin-film data suggests that thin-film experiments are 
indeed relevant probes of the length scale of bulk coopera¬ 
tive dynamics, that may exist independently of any struc¬ 
tural length scale in the material. This approach may be 
rehned with additional cooperative processes, and could 
be adapted to the cases of attractive substrates, free¬ 
standing films, or other geometric confinements, whose 
effects may be crucial on the measured glass-transition 
temperature. 

Acknowledgments 

The authors thank Yu Chai, Olivier Dauchot, Jack 
Douglas, Mark Eton, Florent Krzakala, and James Sharp 
for insightful discussions. This research was supported in 
part by NSERC of Canada and the Perimeter Institute 
for Theoretical Physics. Research at Perimeter Institute 
is supported by the Government of Canada through In¬ 
dustry Canada and by the Province of Ontario through 
the Ministry of Economic Development & Innovation. 


* Corresponding author: thonias.salez@espci.fr 
^ Corresponding author: jforrest@perimeterinstitute.ca 

[1] M. D. Ediger and P. Harrowell, “Perspective: Super¬ 
cooled liquids and glasses,” The Journal of Chemical 
Physics, p. 137, 2012. 

[2] J. H. Gibbs and E. A. DiMarzio, “Nature of the glass 
transition and the glassy state,” The Journal of Chemical 
Physics , vol. 28, p. 373, 1958. 

[3] R. Mari, F. Krzakala, and J. Kurchan, “Jamming ver¬ 

sus glass transitions,” Physical Review Letters, vol. 103, 
p. 025701, 2009. ^ 

[4] H.-N. Lee, K. Paeng, S. F. Swallen, and M. D. Ediger, 
“Direct measurement of molecular mobility in actively 


deformed polymer glasses,” Science, vol. 323, no. 5911, 
pp. 231-234, 2009. 

[5] F. Krzakala and L. Zdeborova, “Glassy dynamics as 
a melting process,” The Journal of Chemical Physics, 
vol. 134, p. 034513, 2011 

[6] E. P. Bernard and W. Krauth, “Two-step melting in 
two dimensions: first-order liquid-hexatic transition,” 
Physical Review Letters, vol. 107, p. 155704, 2011. 

[7] G. Biroli and J.-P. Bouchaud, “Clearing up the mysteries 
of glassy dynamics,” Physics, vol. 6, p. 128, 2013. 

[8] W. Gotze and L. Sjogren, “Clearing up the mysteries of 
glassy dynamics,” Rep. Prog. Phys., vol. 55, p. 241, 1992 

[9] P. W. Anderson, “Through the glass lightly,” Science, 
vol. 267, p. 1615, 1995. 

[10] L. Berthier, G. Biroli, J.-P. Bouchaud, L. Cipeletti, 
and W. van Saarloos, eds.. Dynamical heterogeneities in 
glasses, colloids, and granular media. Oxford University 
Press, 2010. 

[11] A. J. Liu and S. R. Nagel, “The jamming transition 
and the marginally jammed solid,” Annual Review of 
Condensed Matter Physics, vol. 1, 2010. 

[12] G. Paris! and F. Zamponi, “Mean-field theory of hard 
sphere glasses and jamming,” Review of Modern Physics, 
vol. 82, p. 789, 2010. 

[13] L. Berthier and G. Biroli, “Theoretical perspective on 
the glass transition and amorphous materials,” Review 
of Modern Physics, vol. 83, p. 587, 2011. 

[14] M. H. Cohen and G. S. Grest, “Liquid-glass transition, 
a free-volume approach,” Physical Review B, vol. 20, 
p. 1077, 1979. 

[15] A. K. Doolittle, “Studies in newtonian flow. ii. the depen¬ 
dence of the viscosity of liquids on free-space,” Journal 
of Applied Physics, vol. 22, no. 12, p. 1471, 1951. 

[16] H. Vogel, “The law of the relationship between vis¬ 
cosity of liquids and the temperature,” Physikalische 
Zeitschrift , vol. 22, p. 645, 1921. 

[17] G. S. Fulcher, “Analysis of recent measurements of the 
viscosity of glasses,” Journal of the American Ceramic 
Society , vol. 8, no. 6, p. 339, 1925. 

[18] G. Tammann and W. Hesse, “Die Abhangigkeit 
der Viscositat von der Temperatur bie unterkiihlten 
Fliissigkeiten,” Zeitschrift fiir Anorganische und 
Allgemeine Ghemie , vol. 156, p. 245, 1926. 

[19] G. Adam and J. H. Gibbs, “On the temperature de¬ 
pendence of cooperative relaxation properties in glass¬ 
forming liquids,” The Journal of Ghemical Physics, 
vol. 43, p. 139, 196A 

[20] C. Donati, J. F. Douglas, W. Kob, S. J. Plimpton, P. H. 
Poole, and S. G. Glotzer, “Stringlike cooperative motion 
in a supercooled liquid,” Physical Review Letters, vol. 80, 
no. 11, p. 2338, 1998. 

[21] J. D. Stevenson, J. Schmalian, and P. G. Wolynes, 
“The shapes of cooperatively rearranging regions in glass¬ 
forming liquids,” Nature Physics, vol. 2, p. 268, 2006. 

[22] E. Donth, “Gharacteristic length of the glass-transition,” 
Journal of Polymer Science: Part B: Polymer Physics, 
vol. 34, p. 2881, 1996. 

[23] C. L. Jackson and G. B. McKenna, “The melting behav¬ 
ior of organic materials confined in porous solids,” The 
Journal of Ghemical Physics , vol. 93, no. 12, pp. 9002- 
9011, 1990. 

[24] J. L. Keddie, R. A. L. Jones, and R. A. Cory, “Size- 
dependent depression of the glass transition temperature 
in polymer films,” EPL (Europhysics Letters), vol. 27, 



7 


no. 1, p. 59, 1994. 

[25] C. J. Ellison and J. M. Torkelson, “The distribution of 
glass-transition temperatures in nanoscopically confined 
glass formers,” Nature Materials, vol. 2, p. 695, 2003. 

[26] M. Alcoutlabi and G. B. McKenna, “Effects of confine¬ 
ment on material behaviour at the nanometre size scale,” 
Journal of Physics: Condensed Matter , vol. 17, p. R461, 
2005. 

[27] M. D. Ediger and J. A. Forrest, “Dynamics near free 
surfaces and the glass transition in thin polymer films: 
A view to the future,” Macromolecules, vol. 47, p. 471, 
2014. 

[28] F. Varnik, J. Baschnagel, and K. Binder, “Reduction 
of the glass transition temperature in polymer hlms: A 
molecular-dynamics study,” Physical Review E, vol. 65, 
p. 021507, 2002. 

[29] Z. Fakhraai and J. A. Forrest, “Measuring the surface 
dynamics of glassy polymers,” Science, vol. 319, p. 600, 
2008. 

[30] Z. Yang, Y. Fujii, F. K. Lee, C.-H. Lam, and O. K. C. 
Tsui, “Glass transition dynamics and surface layer mobil¬ 
ity in unentangled polystyrene hlms,” Science, vol. 328, 
no. 5986, pp. 1676-1679, 2010. 

[31] Y. Chai, T. Salez, J. D. McGraw, M. Benzaquen, 
K. Dalnoki-Veress, E. Raphael, and J. A. Forrest, “A di¬ 
rect quantitative measure of surface mobility in a glassy 
polymer,” Science , vol. 343, p. 994, 2014. 

[32] H. Yoon and G. B. McKenna, “Substrate effects on glass 
transition and free surface viscoelasticity of ultrathin 
polystyrene hlms,” Macromolecules, 2014. 

[33] K. Ngai, A. Rizos, and D. Plazek, “Reduction of the glass 
temperature of thin freely standing polymer hlms caused 
by the decrease of the coupling parameter in the coupling 
model,” Journal of Non-Crystalline Solids, vol. 235-237, 
p. 435, ism 

[34] D. Long and F. Lequeux, “Heterogeneous dynamics at 
the glass transition in van der waals liquids, in the bulk 
and in thin hlms,” European Physical Journal E, vol. 4, 
p. 371, 2001. 

[35] S. Herminghaus, K. Jacobs, and R. Seemann, “The glass 
transition of thin polymer hlms: some questions, and a 
possible answer,” European Physical Journal E, vol. 5, 
p. 531, 2001. 

[36] J. E. G. Lipson and S. T. Milner, “Percolation model 
of interfacial effects in polymeric glasses,” European 
Physical Journal B, vol. 72, p. 133, 2009. 

[37] S. Mirigian and K. S. Schweizer, “Communication: Slow 
relaxation, spatial mobility gradients, and vitrihcation 
in conhned hlms,” The Journal of Chemical Physics, 
vol. 141, p. 161103, 2014. 

[38] J. Baschnagel and F. Varnik, “Computer simulations of 
supercooled polymer melts in the bulk and in conhned ge¬ 
ometry,” Journal of Physics: Condensed Matter, vol. 17, 
no. 32, p. R851, 2005. 

[39] B. A. Pazmino Betancourt, J. F. Douglas, and F. W. 
Starr, “String model for the dynamics of glass-forming 
liquids,” The Journal of Chemical Physics, vol. 140, 
p. 204509, 2014. 

[40] Z. Zhang, P. J. Yunker, P. Habdas, and A. G. Yodh, “Co¬ 
operative Rearrangement Regions and Dynamical Het¬ 
erogeneities in Colloidal Glasses with Attractive Ver¬ 
sus Repulsive Interactions,” Physical Review Letters, 
vol. 107, p. 208303, 2011. 

[41] J. Mattsson, H. M. Wyss, A. Fernandez-Nieves, 


K. Miyazaki, Z. Hu, D. R. Reichman, and D. Weitz, “Soft 
colloids make strong glasses,” Nature , vol. 462, p. 83, 
2009. 

[42] A. Donev, 1. Cisse, D. Sachs, E. A. Variano, F. H. Still- 
inger, R. Connelly, S. Torquato, and P. M. Chaikin, “Im¬ 
proving the density of jammed disordered packings using 
ellipsoids,” Science , vol. 303, no. 5660, p. 990, 2004. 

[43] P. Pal, C. S. O’Hern, J. Blawzdziewicz, E. R. Dufresne, 
and R. Stinchcombe, “Minimal model for kinetic arrest,” 
Physical Review E, vol. 78, p. 011111, 2008. 

[44] A. S. Keys, A. R. Abate, S. C. Glotzer, and D. J. Durian, 
“Measurement of growing dynamical length scales and 
prediction of the jamming transition in a granular mate¬ 
rial,” Nature Physics , vol. 3, p. 260, 2007. 

[45] S. L. Simon, J. W. Sobieski, and D. J. Plazek, “Vol¬ 
ume and enthalpy recovery of a-polystyrene,” Polymer, 
vol. 42, p. 2555, 2001. 

[46] A. Sahnoune, F. Massines, and L. Pich, “Ultrasonic mea¬ 
surement of relaxation behavior in polystyrene,” Journal 
of Polymer Science, Part B: Polymer Physics, vol. 34, 
p. 341, 1996. 

[47] S. Kahle, J. Korus, E. Hempel, R. Unger, S. Horing, 
K. Schrotter, and E. Donth, “Glass-transition coopera- 
tivity onset in a series of random copolymerspoly(n-butyl 
methacrylate-stat-styrene),” Macromolecules, vol. 30, 
p. 7214, 1997. 

[48] M. L. Williams, R. F. Landel, and J. D. Ferry, “The tem¬ 
perature dependence of relaxation mechanisms in amor¬ 
phous polymers and other glass-forming liquids,” Journal 
of the American Chemical Society, vol. 77, p. 3701, 1955. 

[49] J. Langer, “Dynamics and thermodynamics of the glass 
transition,” Physical Review E, vol. 73, no. 4, p. 041504, 
2006. 

[50] C. A. Angell, “Formation of glasses from liquids and 
biopolymers,” Science , vol. 267, p. 1924, 1995. 

[51] J. Guo, X. Bian, Y. Zhao, S. Zhang, T. Li, and C. Wang, 
“Correlation between the fragility of supercooled liquids 
and thermal expansion in the glassy state for gd-based 
glass-forming alloys,” Journal of Physics: Condensed 
Matter, vol. 19, p. 116103, 2007. 

[52] S. F. Edwards and T. Vilgis, “The dynamics of the glass 
transition,” Physica Scripta, vol. T13, p. 7, 1986. 

[53] J. Zhao, S. L. Simon, and G. B. McKenna, “Us¬ 
ing 20-million-year-old amber to test the super- 
arrhenius behaviour of glass-forming systems,” Nature 
Communications, vol. 4, p. 1783, 2013. 

[54] D. M. Coined, G. B. McKenna, F. J. J., A. Lee, D. B. 
Curliss, K. B. Bowman, and J. D. Russel, “Isochoric and 
isobaric glass formation: Similarities and differences,” 
Journal of Polymer Science Part B: Polymer Physics, 
vol. 35, p. 1561, 1997. 

[55] J. C. Dyre, “Colloquium: The glass transition and elas¬ 
tic models of glass-forming liquids,” Reviews of Modern 
Physics , vol. 78, no. 3, p. 953, 2006. 

[56] F. Lechenault, O. Dauchot, G. Biroli, and J. P. 
Bouchaud, “Critical scaling and heterogeneous superdif¬ 
fusion across the jamming/rigidity transition of a granu¬ 
lar glass,” Europhysics Letters, vol. 83, no. 4, p. 46003, 
2008. 

[57] L. Berthier, “Finite-size scaling analysis of the glass tran¬ 
sition,” Physical Review Letters, vol. 91, p. 055701, Aug 
2003. 

[58] C. Donati, S. C. Glotzer, P. H. Poole, W. Kob, and S. J. 
Plimpton, “Spatial correlations of mobility and immo- 



bility in a glass-forming Lennard-Jones liquid,” Physical 
Review E , vol. 60, no. 3, p. 3107, 1999. 

[59] I. Karatzas and S. E. Shreve, Brownian motion and 
stochastic calculus. Springer, 1988. 

[60] J. A. Forrest and K. Dalnoki-Veress, “When does a glass 
transition temperature not signify a glass transition?,” 
ACS Macro Letters , vol. 3, no. 4, pp. 310-314, 2014. 

[61] C. B. Roth and J. R. Butcher, “Glass transition 
and chain mobility in thin polymer films,” Journal of 
Electroanalytical Chemistry, vol. 584, no. 1, pp. 13 - 22, 
2005. 

[62] A. N. Raegen, M. V. Massa, J. A. Forrest, and 
K. Dalnoki-Veress, “Effect of atmosphere on reduc¬ 
tions in the glass transition of thin polystyrene films,” 
European Physical Journal E, vol. 27, p. 375, 2008. 


[63] J. Brandrup, E. H. Immergut, and E. A. Grulke, Polymer 
Handbook . Wiley-Interscience, 1999. 

[64] M. Rubinstein and R. H. Colby, Polymer Physics. Oxford 
University Press, 2003. 

[65] P. Scheidler, W. Kob, and K. Binder, “The relaxation 
dynamics of a simple glass former confined in a pore,” 
EPL (Europhysics Letters) , vol. 52, no. 3, p. 277, 2000. 

[66] B. Frick, R. Zorn, and H. Buttner, eds., International 
Workshop on Dynamics in Confinement. J. Phys. IV 10, 
Pr7, 2000. 

[67] J. A. Forrest, “What can we learn about a dynami¬ 
cal length scale in glasses from measurements of surface 
mobility?,” The Journal of Chemical Physics, vol. 139, 
p. 084702, 20 It 



